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1. Introduction 

The progress of the string program as a theory of quantum gravity and other interactions 
is currently impeded by a fundamental difficulty, namely we do not know how to formulate 
string theory in a time-dependent background in general, and how to understand many 
issues related to cosmology such as the origin of our universe and the nature of dark energy 
in particular. This baffling situation leads to lots of debates about whether string theory 
is the correct theory of nature, and whether string theory has any predictive power if there 
exists a vast landscape of meta-stable vacua. It goes without saying that string theory 
has been tremendously successful in resolving some of deeper conceptual problems such 
as whether gravity is compatible with quantum mechanics, but only in some unrealistic 
backgrounds such as a flat background and an anti-de Sitter background. In some cases, 
we even have a non-perturbative formulation, for instance, a CFT is a non-perturbative 
theory in the AdS/CFT duality. Nevertheless, until we have a theory for time-evolving 
backgrounds, string theory can not claim to be the theory of our universe. 

We shall not try to attack the ultimately difficult problem of formulating string theory 
in a general or even an "on-shell" time-dependent background in this note. Our purpose 
is rather pragmatic, we will try to work out part of string quantization in a de Sitter 
background, with applications to inflation as well as to a later universe dominated by dark 
energy in mind. For instance, we would like to know how different the spectrum of a string 
in the de Sitter space is from that in the flat spacetime, whether this spectrum enables 
string production during inflation. If not, whether strings are created at the end of inflation 
when the Hubble constant undergoes transition from a constant to a decreasing function. 
In the later universe such as the current epoch, our universe is again dominated by energy 
of almost constant density. Although a cosmic string, if exists, is largely governed by 
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classic dynamics, it is certainly of interest to know whether its spectrum is modified in 
some extremal limit. 

The answers to the above questions seem to be yes. As we shall see, the dynamic 
equation for the field corresponding to a fixed state contains a new term induced by "string 
mode creation" (to be explained shortly). This term depends on the Hubble constant, thus 
it renders string creation possible in the end of inflation. This term begins to be comparable 
to the usual term in the "mass" spectrum when the oscillation numbers are large enough. 
We put mass into quotation marks since there is no notion of mass in a de Sitter space. 

Note that, we will exclusively deal with "small" string states in this paper, by a small 
string we mean that the string modes are mostly oscillating. It is known that cosmic 
strings are "long" strings, namely the dominating modes are not oscillating modes so that 
the major part of the string co-moves with the expansion of the universe. 

There is a series of papers on first quantizing string in de Sitter space by de Vega 
and Sanchez and their collaborators, see for instance (3j-[f|. In their work, they fix all 
the degrees of freedom of the world sheet metric, and find an exact classical solution. To 
quantize, they propose two methods, one is to quantize the fluctuation around an exact 
solution of the center of mass JlJ], and the other is to propose a quantization condition 
semiclassically ||. 

In this note, we propose a new approach. We leave one degree of the world sheet 
metric unfixed, and then eliminate it by a constraint. To the leading order, our result 
agrees with that of de Vega and Sanchez. However, there is a subtle difference: our method 
is approximate in choosing the gauge. Despite of this, the answer is exact once the gauge 
is chosen. The most important consequence of our result is that the mode creation and 
annihilation operators are still time-dependent even after diagonalization. Thus, a state 
created by these operators is itself time-dependent, and this will have some effects in the 
dynamic equation for the corresponding field. 

We will present our approach in the next section and carry out the first quantization 
in Sect. 3. We will discuss possible applications to inflation and obtain some conclusions in 
the last section. A discussion on the mode creation on a string is left in the appendix. 

2. Action and gauge choice 

We start with the Polyakov action in a general background 

S = "I^ 7 / ^^^^^d^G^iX) , (2.1) 

where h = deth a f, (a, b run over values (r, cr).), < a < 2n, and G^ u is the string frame 
metric. Here we suppose that the dilaton is constant, and G^ (/i,^ run over 0, 1, 2, 3. ) 
is the metric of de Sitter space in comoving frame 

ds 2 = -dt 2 + e 2Ht (dx*) 2 , i = 1, 2, 3 . (2.2) 



Of course, if we naively take ( [2.1| ) as the whole story, then in the background (2^) a 



quantum string is not well-defined, since we do not have a two dimensional conformal field 
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theory. We shall assume that there is a hidden sector making the whole world-sheet action 
conformally invariant (as an example, in a merely AdS space, the bosonic string action is 
not conformally invariant, and we need another sector from a sphere as well some other 
terms due to flux and fermionic degrees. As for a de Sitter space, we assume a scenario 
such as KKLT compactification making the story complete). 

By means of the classical world-sheet symmetry, we can set determinant of the world- 
sheet metric to —1 after choosing the temporal gauge and diagonalizing world-sheet metric. 
Following this strategy, there is only one component of the world sheet metric h aa left 
unfixed, which is non-dynamical. We assume that h aa depends only on time, then the 
target space coordinates X 1 can be solved in terms of h aa . When we quantize the field 
X 1 , h aa will be promoted to an operator. In order to obtain the on-shell condition, we 
then impose the constraint from the variation of h aa on physical states. In more detail, 
the constraint is the integral of the variation of h aa due to our assumption. 

Now we perform the steps summarized above. By choosing a proper gauge, we can fix 
the redundancies in the Polyakov action, and make the equations of motion simple. Set 

T = t , h™ = , -h = 1 . (2.3) 
where t is the comoving time. Under this gauge choice, the action becomes 

There are two independent constraints due to functional variation of h ab , that is, 



1 - e 2Ht d t X\t, a) d„X% a) = , (2.5) 



Aira 



1 -{-feLMj + e^fcMjfftl'M)) 2 !^!'^)) 2 ]} = . (2.6) 



87TO 

The equations of motion corresponding to the functional variation of X^ are 

d t (e 2Ht h CT(T (t,a)d t X i (t,a))-e 2Ht d (T (h-^t,a)d (T X i (t,a))=0 , (2.7) 
dtK^a) + H e 2Ht [h aa (t,a)(d t XHt,a)) 2 - h-^a^X'^a)) 2 } =0 . (2.8) 



The second equation of motion ( p.8[ ) comes from the functional variation of X°, which 
is non-dynamical according to our gauge choice. This becomes another constraint, but for- 
tunately it can be derived from the other three equations. So at last we get two constraints 
and three equations of motion (each for one spatial coordinate) . 

The conjugate momentum of X l (t, a) is given by IT(i, a) = e 2Ht h aa (t, a) dtX l (t, a) 
The Hamiltonian is then 

p 2Ht p u r± \ 

E = W J da[ e^Ht +haa(t,a)(d t X%a)) 2 + h^(t,a) (d a X^t,a)) 2 } (2.9) 

~ 2^7 / dah ™( t > <T ) ■ ( 2 - 10 ) 
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In the last step we used (|2.6|) , which is satisfied only by physical states. We use the 
symbol ~ instead of = to show that the equality is satisfied only by physical states. Hence 
J dah a(J is in fact just the energy of a physical state with respect to comoving time. 
For simplicity, we will set a' = 1 hereafter. h aa is non-dynamical, and is determined by 
X l (t,a) through ( p.6[ ). If we eliminate h aa , (2.7) will become nonlinear equations which 
is hard to solve. Instead, we will first treat h aa as an independent variable, solve the 
equations of motion for X l (t, a) in terms of h aa and then fix h aa by ( |2.6| ). In order to solve 
the equation of motion (|2.7| ), we make an assumption that 

h aa (t,a)=u(t). (2.11) 

This is the only approximation we execute in this note. For strings oscillating fast, tu{t) 
may be viewed as an average of h aa (a,t) along a. 

Hereafter we will just write lo instead of Lo(t) for simplicity, but keep in mind that u> is 
in fact a function of time. Also, we caution that upon quantization, X 1 become operators, 
so does to, namely u> is not to be viewed as a usual function. 

For physical states, E ~ u is the energy of the string, the equations of motion become 

d^ 2 d t X% a)) - uj- 2 7]- 2 d 2 a X l (t, a) = , (2.12) 
where r] = ■ A general solution is 

X l {t,a) =X + f durj 2 {u) P i + 7?(t)[ =e -iJ t du\ m (u) e ima 

J meZ/{0} V 2 \ X m{t)\ 

I d m (t) ft du\ m (u) g-imcn ^ \£\ 

V2\KAtT\ 

al(t) =^f~m{t)e 2 ^ td ^\ alit) = ^/-mit) e^/'^W .(2.14) 



where we have definedA m = sgn{m)\J — ?]9 t 2 (r/ _1 ) = sgn(m)y — (H + ^j) 2 — <9t(^j), 
where the function sgn{m) = 1 for m > 0, and sgn(m) = — 1 for m < . We will work in 
the situation that A m s remain real, which means that the string is oscillating in time. This 
is to be dubbed as a small string, since it is not stretched too much with the expansion of 
the universe. 

The real condition for X m is that - (H + §^) 2 - d t {^) > 0. The Hermiticy of 
X l (t, a) requires that {a l m )^ = o!_ m and (a l m )^ = a l _ m . Thus the conjugate momentum of 
X 1 becomes 

rr^a) =—{ P i + Y <M ^^1^ 

| d m (t) ft du\ m (u) e -im<Ti^ 

To quantize, impose the equal time canonical commutation relations 

[X i {t,v),Xi{t,o')] = [U t (t,a),U j (t,a')}=0 , (2.15) 
[X i {t,o),Il j {t,a')}=i5 j 5{<j -a') . (2.16) 
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which are equivalent to imposing the commutation relations 



[x\ x j ] = \p\p>] = , [x\p>] = i5 ij (2.17) 

771 ■ ■ 

[a l m ,a 3 n ] = , [a l m ,a J n ] = —5 l:, 5 m ~n • (2.18) 



Note that a l m and a % m depend on u implicitly. Thus when we impose (|2.1q ), we have 
also promote u to be an operator, which commutes with other operators including to. To 
simplify the notations, we will not distinguish operators and functions explicitly except 
that u, p and N are assumed to be operators when constructing states. 

3. Quantization 

We now study constraints (2.5), (EO|) and (|2.8|). fl2.8|) can be derived from (|2.5|), (p.6|) 



together with the equation of motion (|2.7|). Thus only ( |2.5| ) and ( |2.6| ) need to be considered. 
According to our assumption h aa (t, a) = u>(t), the most important parts of these constraints 
are their average over a. Thus the constrains become, 



P = j ^e 2Hi ■■d t X%&)d t ,X\t,&) 



2. ^ 2uj { a ~ m am ^-m^m} — ! (3-1) 



m>0 



and 



* = J ^{hL[-l + e 2Ht (dtX*(t,a)) 2 ]+e 2Ht (d a X\t,<j)) 2 } (3.2) 
_ ^ | (P 1 ) 2 | y tj 77 2 2 m 2 , < , , 

m>0,i 

+ r [( 3 + , Am)2 + ^ a ,_ m ~ alm ^ (3 3) 

There are also infinitely many constraints corresponding to positve modes of ( |2.5| ), ( |2.C| ) 
in the Fourier expansion in terms of a. Note that our temporal gauge choice only eliminate 
the temporal degrees of freedom. These infinitely many constraints, whose analog in flat 
spacetime are the conditions L n = L n = 0, will further eliminate the unphysical degrees of 
freedom due to the longitudinal excitations. In this note we will not discuss them in detail, 
and will just focus on the mass shell condition, whose analog in flat spacetime is Lq = 0. 
One may ask that why there are still two set of constraints while only the longitudinal 
degrees of freedom need to be eliminated. Are our constraints too strong? The answer 
is no, because we have used the assumption that h aa is independent of a, which imposes 
another constraint. If we re-introduce h aa as an arbitrary function of <r, the counting of 
degrees of freedom will be correct. 

Define the occupation number n l m = a i _ m d l m , and the level n = Ylmi mn mi an< ^ 
similarly for n. Then the vanishing of ( |3.1| ) on physical states is just the level matching 
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condition n = h. Note that the condition ( |3,l[) is also the translational invariance condi- 
tion along a. We can make a linear transformation to define another set of creation and 
annihilation operators, 

\a m \ 2 - \P m \ 2 = \a m \ 2 - m \ 2 = 1 , (3.5) 
a m f3- m = Oim$-m • (3.6) 

The conditions ([II]) and ([O]) ensure that the new operators satisfy the same commutation 
relation as ( 2.1 8| ) . The above transformaion may be viewed as Bogoliubov transformation 
on the world-sheet. The most general form is 

a m = cosh( 7m ) e iS ™ +i4> ™,a m = cosh( 7m ) e iSM ™ , (3.7) 
P m = sinh( 7m ) e 1 ^ , p m = sinh( 7m ) e**" . (3.8) 



where r y m ,<p m ;''Pm an d are real. Constraint ( p.l| ) remains the level matching con- 
dition, but now in terms of the new occupation number operator N = N, defined by 
N = ^2 im mNm, and = A l _ m A l m , and similarly for N. 

We now choose the parameters 7m and 5 m properly to diagonalize the constraint (|3,2|). 
The conditions are 

2 

cosh 2 ( 7m ) + sinh 2 ( 7m ) = — ^— [(^) 2 + \ 2 m + ^] , (3.9) 

sinh(2 7m ) = [(JL - i\ m ) 2 + e-^S'duXm^-iSm (3 1Q) 

2m A m 7/ w z 

where 5 m is chosen to make sinh(2 7m ) real, and this can be done. 
In terms of the new creation and annihilation operators 

m>0,i 

- 4 + Ae 2Ht + + 2 - ' 4 • 

The fact that we need to form a Bogoliubov transformation implies that if we start with 
a state constructed by the original operators a l m , there will be mode creation along the 
string at a later time. 

We now discuss implication for dynamics of fields viewed as coefficients in the expansion 
of a general string state. 

Since d% is a Killing vector, p' 1 is conserved, and we will just work in momentum 
representation. Any physical state should satisfy the constraint 34? = 0, and can be 
expanded in terms of common eigenstates of the occupation number operator N^s, N^s, 
p l as well as u. That is, 

\4>>= Yl WlNl, Nl,N l m , X^iN^N^uJ) , (3.12) 
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where the eigenvalues labeling the states must satisfy the relation —u 2 +^^t+2(2N+Eo) = 
0, and N = N. 

When writing down the action for the above general state, it is important to keep 
in mind that the inner product involves an integral over space so there is a nontrivial 
Hermticity condition. For instance, we consider a scalar particle with mass to, whose wave 
function must satisfy 

(□ - m 2 )<t){x) = \-d 2 t - 3Hd t + e~ 2Ht (d*) 2 - m 2 ]0(x) = . (3.13) 

Written in momentum representation, the last two terms are —e~ 2Ht (p 1 ) 2 — m 2 = —E 2 , 
where E is the comoving energy. 

In string field theory, we assume that the action of string state has the form of 

S = j dt< (f>\d? + 3Hd t + E 2 + AJT|0 > (3.14) 

where A is the Lagrangian multiplier. In writing down this action, we have only consid- 
ered the mass shell condition and have omitted other constraints corresponding to the 
positive Fourier modes. To treat the problem more rigorously, one should introduce more 
Lagrangian multipliers. From this action, we can easily get the evolution equation of the 
string state. By variation of \<j> >, we have 

(d t 2 + 3H8 t + E 2 )\cp >= [8 2 t + 3Hd t + e~ 2Ht (p*) 2 + 2(N + N + £ )]|</> >= . (3.15) 

This equation contains an unphysical component which is to be discarded due to the 
fact that in the action the inner production automatically projects out the unphysical 
component by imposing the constraint Jt?\(f) >= 0. 

We explain some subtleties in deriving this second order equation. Because of the non- 
flat metric, the measure of the integral volume is dx 3 \/—G = dx 3 e 3Ht , so the inner product 
should be defined as J d 3 xe 3Ht (j){x)*(f)(x). With this inner product, d 2 is not Hermitian. To 
get a Hermitian operator, we should replace d 2 with df + 3Hdt- There is no addition term 
caused by polarization indices since all the creation operators are properly normalized. 

Now, different excitation modes of the string correspond to different particles in space- 
time, thus the coefficient </>(iV m , iV m , cj,p l ) is the wave function of single particle. From 

now on, we will use notation |iV m , N^u^p 1 > instead of \Nl,Nl, _/V m ,iV m , io,p l > 

for simplicity. Since the basis |iV m , N^^uj^p 1 > evolves with time, this dependence on time 
will be transmitted into the equation of motion for <fi through the action of d 2 + 3Hdt ■ 



Using the differential equation for a l m and a m , (|2.14[) , and the definition of A l m , A l m , we 
have 

41 _ „ 4« _i_ J Al A _ 7. At _i_ A} (IIP,) 

—m m ' ""m -^—mi —m °m m ' "wi ■ rL — m ' 

r _ J{<t>-i>) S r* a* _ ft* A* i \ n * 2 p -2i f du\ m (u) _ o*2 2i f du\ m (uh_^rn_\ 

a j£ • 2m f tt _i u> \ 

e (-i5-i4>-ij>) jj Ut 2uj 1 ui v-" 2ui> ^ -j^y-j 

\^)2(H+^ + [d t (%)] 2 ' 
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c m = e^-^{a*J* m - fa&l + [(d*) 2 e" 2i J" duX ^ - (^V^* 1 ^)]-^} 

a _uj 2m / tt _i Oj_ \ 

e (-iS-i<f>-iil>) ut 2uj 1 u> v 17 2ui> ^ ^g-j 

J —J{4>-^){^ a* - R R* I ^ m U /?* p 2i f du X m (u) _ * a -2i f du A m (u)ii 

H + AHm , 9 , r « w . cj . „ cj 4m 2 w , , 



" 2AS l (^-^-M Am) l w m mVt 2u 2u ,l 2u 2u 

-i<j) , (3.19) 

d m = e^\a m a m - fr$ m + ^[a m p m e 2i ? duX ™^ - a m (3 m e~ 2i f 

H -\- TT— 4Hm , 9 , , „ 9 Co r cj . . Co Am 2 Co , , 

{ A^ + A^d 2 — + 2(tf + — — + ^— ]} 



" 2A^-a^-^A m ) 1 - ™' "^2. ^ 2,' ^ ufilu, 
—zip . 

Note that |c m | = if. We can choose and Y> to set d m = d m = 0, and then the phase of 
c rn is also fixed. 

The eigenstates of and N m are just \\ m ^(A i _ m ) N ™(A i _ rn ) N ™ \Q,u >, where \Q,w > 
is defined as A m \Q,io >= 0, A m \Q,L0 >= for all m,i. From \Cl,co > +A m d t \Q,io >= 
0, we get 

d t \n,u>=-^<&AL m AL m \il,u,> , (3.20) 

and 

a t |AT^iV^ W >=^[ Cm ^Ij n -C^: m lLJ|iV n fc ,iV n fe ,u;> , (3.21) 

d 2 t \Nt N^co >=-J2 (! + 2iV ™ + < + ^K*' ^ > . ( 3 - 22 ) 
+{— ^ 2 c m c* t A l _ t A l _ t A J m A J m + ^2 c m A m A l m — ^ c* m A % _ m A l _ m (3.23) 

m,i,l,j m,i m,i 

+ c ™ A\ A\ Ai Ai + c\ c* m AU AU AL m Ai ro ]}|JV* N k , to > , (3.24) 

m,i,l,j 

tt 2 -2Ht 

Co\NlN k n ,uo>= V - > , (3.25) 

CO 

N + ti = Y, 2m(c m A m A m + c* m A\ m ~A\ m ) . (3.26) 

m,i 

From the explicit expression above, we see that the off diagonal parts of dt\N, N,uj > and 
®t I N m , N m , to > are unphysical if \N m , N m ,co > is physical. So the off diagonal part is 
orthogonal to physical states. 

Just considering the physical part of the following equation 

{3 t 2 + 3Hd t + (p*) 2 e- 2Ht + 4N + 2E O }\0 >= , 
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we have 



{8 2 t + 3Hd t + (p 1 ) 2 e~ 2Ht + AN + 2E Q 

-^H 2 (l + 2< Ni n + N l m + iV^)}0(iV n fc , NH,p) = 0. (3.27) 

m,i 

The above equation is the main result of this note. We see that in addition to the term 
4iV, there is an additional term which is quadratic in creation numbers with a prefactor 
H 2 . This term could become comparable to the linear term N for a fixed H 2 (measured 
in the string unit since we have set a' = 1). Restoring the string scale M s , we find that 
the new term is comparable to the old term M 2 N when N ~ M 2 /H 2 . This is of course a 
large number during inflation because H is much smaller than the string scale, if we hope 
that an effective field theory is valid. 

This quadratic term is negative when viewed as a contribution to the mass squared 
m 2 , thus it appears possible to have an effective negative mass squared if the quadratic 
term becomes dominating. This will never become possible, since in our approach so far 
we have assumed real X m in the mode expansion and we are dealing with "small" 

strings mostly oscillating in time. It can be checked that as long as the real condition 
on A m is met, the quadratic term in (|3.27| ) will never make the mass squared negative. 
Nevertheless, the fact that this term reduces the mass squared comes as a surprise. We 
may just imagine that this is a quantum effect for a highly excited state due to the mode 
creation on the string. We leave a discussion on mode creation to the appendix. 



4. Comparison with earlier results 

The spectrum we get is 

a'M 2 = 4N + 2E -Y J H " i 1 + 2N L < + K. + Ki), t 4 - 1 ) 



which can be read from ( 3.27 ). The previous result obtained by the method of (TJ is 



/ „ ,2 , V- 2n 2 - H 2 M 2 a' 2 Ar 2 - H 2 M 2 a' 2 , A , 

a'M 2 = 24 > + 9 AT , 4.2) 

^ Vn 2 - H 2 M 2 a' 2 y/l - H 2 M 2 a' 2 

where we cite the formula in the form appearing in ||, and the dimensionality of the 
spacetime is 25 there. We can see that when H = 0, both results return to the flat 
spacetime spectrum, 4iV + zero point energy. The difference between the zero point energy 
is due to the different dimensionality. When H ^ 0, both results have a term proportional 
to H 2 , but the coefficients are different. What's more, our result shows that the spectrum 
depends not only on the level N, but also on the specific excitation. This difference 
may due to the gauge choice of the two approaches. As we have mentioned before, we 
fix the worldsheet time r, set the determinant of the worldsheet metric to be one, set 
one component of the worldsheet metric h 12 = 0, and leave another worldsheet metric 
component h aa unfixed. To set the determinant of the worldsheet metric to be one, we 
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have used the classical conformal symmetry of the action, which does not exist in the whole 
theory. While in the earlier approach, for example, [jlj], they choose the conformal gauge and 
leave the worldsheet coordinate r unfixed. Since there is no conformal symmetry, there is 
no guarantee that we should get the same result. The other approach previously developed 
in namely, the semiclassical quantization, also take the same conformal gauge, and make 
a circular string ansatz, and then they propose a quantization condition. The spectrum is 
5.9n,n G Nq in [Q]. In 0|, the authors compare this result with that obtained in 



/™2 



a m 



|1[ by calculating the maximum excitation number of string states in de Sitter spacetime. 
In m, the maximum number of a single excitation is N max = jjt^t, which obtained in 
|H is N max = . In our approach, the condition for the state to oscillate is > 0. 
Together with ( 3.12Q , we will get roughly N max = for all the oscillating modes to have 
real frequency. Our result is slightly larger than the previous results. One possibility is 
that our solution is more general than the previous ones. For instance, in [||, only circular 
solutions are considered, and in only expansion around an exact solution is considered. 
While here in our paper, only one approximation is made, namely, the worldsheet metric 
h aa does not depend on a. Thus we might have found more solutions. 



5. Discussions and conclusion 

In this paper, we first quantized a general oscillating string in a de Sitter space , with the 
only approximation that h au depends only on time. This quantity becomes the energy 
density along the string after we impose constraint condition on it. So our approximation 
amounts to averaging the energy density along the string. Aside from this, our treatment 
is exact. 

Apparently, our main result ( 3.27] ) differs from the old result by a negative contribution 



to the mass squared of the string. This term which is quadratic in oscillation numbers 
increases quickly when we consider more highly excited states. Applying this result to 
inflation, there is virtually no physical effect during inflation except for modification of the 
mass spectrum, for this new term depends only on the Hubble constant thus remains a 
constant for a given state. It is a well-known result that for a particle of constant mass, 
there is no particle production during inflation. Although we have worked with a constant 
Hubble constant, some of our results can be generalized to a non-constant Hubble constant. 
We expect string creation in the end of inflation can be induced by this new term , for the 
Hubble constant is no longer a constant in this short reheating period . A similar effect 
is discussed in Jjj and || except for ad hoc coupling to some moduli. We leave a detailed 
investigation of string creation to a future work §. 

In this paper we have restricted attention to "small" string states, namely strings with 
oscillating modes only. It can be expected that the phenomenon of string mode creation 
and its induced effects on the equation of motion prevails for the long strings which stretch 
with the expansion of the universe. Again, detailed result will be presented in §. 

We also expect the new term we discovered will have some effects for cosmic strings 
at later times. 
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6. Appendix 

To investigate the question of mode creation, we need to consider the expectation value 
of the occupation number operator on a state that does not change with time, that is, 

the eigenstates of the time independent number operator. Because A± m and A± m do not 

vanish, the corresponding number operators N and N depend on time. Thus we need to 

find a set of creation and annihilation operators b l m and b l m , s.t. b % ± m = b± m = 0. b± m and 
b l ± m are linear combinations of A± m and A± m . As we have mentioned under fl3.4|) , the 
most general form of linear transformation preserving the commutation relation is to set 

V m = cosh( 7 ;je^ + <^ + sinh(yje^I*_ m , (6.1) 
V m = cosh( 7 ;je^ + ^Ij n + sinh( 7 ^)e^^_ m . (6.2) 

Demanding that 

bin = W'm + O cosh(yje^+^ + j' m sinh( 7 ;je^+<* + c m sinh( 7 Je^]4{6.3) 
+ [i# nS inh(y m ) e i #» +y TO cosh(yje i ^ +Ccosh( 7 ;je^ + ^]i j : m = , 

C = IW'm + D cosh(yj e ^+^ + i m sinh(yje^+^ + c m sinh( 7 Je^H 1^6.4) 
+ sinh( 7 ^)e^ + j' m cosh( 7 ^)e^ + c* m cosh^'Je^+^A^ = . 

Then we have 

j' m + M(c m e- iS '™) = , (6.5) 
5' m + [tanh( 7 ^) + coth(yj]9(c m e- i5 ™) = . (6.6) 

Suppose that at time t = to, the two sets of operators are identical, e.g. b l m = A l m (tQ), b\ n = 
A l m (to). Then the Hilbert space spanned by eigenstates of N' = 6!_ m &m and N' = &L m &m 
is identical with that of N(to) and N(to). Thus we have the initial condition = 8' m = 0. 
Denote \t Q > a state that does not change with time, then < to\N(t)\to > represents the 
change of the total mode number. 

<t \N + N\t >= -^2mHsmh(2"f' m )cos(5' m )(N l m + N l m + l) . (6.7) 

m,i 

When p = 0, both the real part and the imaginary part of c m is pure oscillating, so the 
average of j' m and 5' m vanish, with the initial condition, we will have n / m ~ 5' m ~ 0. Then 
< t \N\t >»0. 
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